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I.  INTRODUCTION 

This  report  summarizes  progress  made  during  the  period  September  1,  1982 
to  August  31,  1983,  under  the  project  entitled  "Ignition  of  Solid  Propellants 
and  Propagation  of  Burning  Propellant  Cracks"  (Contract  No.  N00014-79-C-0762) . 

The  overall  objective  of  this  investigation  is  to  achieve  better  under¬ 
standing  of  the  ignition  mechanism  under  transient  pressure  loading  conditions 
and  to  investigate  propagation  of  burning  solid  propellant  cracks.  Specific 
objectives  of  this  study  are: 

1.  to  observe  the  ignition  phenomena  of  nitramine  composite  propellants,  and 
to  measure  ignition  delay  under  various  operating  conditions; 

2.  to  compare  ignition  delays  of  different  types  of  propellant  samples  under 
identical  flow  conditions; 

3.  to  study  the  effect  of  bindef  and  oxidizer  ingredients  on  ignition  by 
comparing  the  ignition  processes  of  propellants  with  different  ingre¬ 
dients; 

4.  to  propose  a  theoretical  model  for  the  ignition  of  nitramine  propellants, 
and  to  gain  a  deeper  understanding  of  the  ignition  process; 

5.  to  observe  crack  propagation  phenomena  in  burning  solid  propellant  sam¬ 
ples; 

* 

6.  to  study  the  mechanism  of  crack  propagation  and  branching; 
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7.  to  determine  the  effect  of  pressurization  rate  and  sample  geometry  on 
crack  propagation  velocity;  and 

8.  to  evaluate  the  ratio  of  effective  mass  burning  rate  to  the  steady-state 
mass  burning  rate  of  the  sample  under  various  operating  conditions. 

The  following  is  a  list  of  papers  published  during  the  past  year  under 

the  support  of  this  contract. 

1.  "Ignition  of  Composite  Propellants  in  a  Stagnation  Region  Under  Rapid 
Pressure  Loading ,"  Nineteenth  International  Symposium  on  Combustion,  Hie 
Combustion  Institute,  1982,  pp.  757-767  by  M.  Kunar,  J.  E.  Wills,  A.  K. 
Kulkami,  and  K.  K.  Kuo. 

2.  "Ignition  of  Nitramine  Propellants  Under  Rapid  Pressurization," 
AIAA/SAE/ASME  19th  Joint  Propulsion  Conference,  AIAA  Paper  83-1194,  Seat¬ 
tle,  Washington,  June  27-29,  1983  by  S.  YU,  W.  H.  Hsieh  and  K.  K.  Kuo. 

3.  "Mechanism  of  Crack  Propagation  in  Burning  Solid  Propellants,"  1983 

JANNAF  Propulsion  Systems  Hazards  Meeting,  Los  Alamos,  tw,  July  13-15, 
1983  by  K.  K.  Kuo,  J.  G.  Siefert,  M.  Kumar,  W.  H.  Hsieh,  and  E.  Andi- 
roglu.  -  - -  -  - - - 

4.  "Crack  Propagation  in  Burning  Solid  Propellants,"  Hie  Proceedings  of  the 
Ninth  International  Colloquium  on  the  Dynamics  of  Explosions  and  Reactive 
Systems,  Poitiers,  France,  July  3-8,  1983  by  J.  G.  Siefert  and  K.  K.  Kuo. 

5.  "Nitramine  Propellant  Ignition  Under  Rapid  Gas  Loading  Conditions,"  20th 
JANNAF  Combustion  Meeting,  Monterey,  California,  Oct.  17-20,  1983  by  S. 
Yu,  W.  H.  Hsieh,  and  K.  K.  Kuo. 
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Besides  the  above  publications,  research  results  obtained  from  this 
project  were  also  presented  at  the  following  workshops. 

1.  "Ignition  of  Nitramine-LOVA  Propellants,"  presented  at  the  JANNAF  Work¬ 
shop  on  Characteristics  of  Ni t r ami ne-LOVA  Gun  Propellants  held  at  NOS, 
Indian  Head,  M3,  June  14-15,  1983  by  K.  K.  Kuo,  S.  Yu,  and  W.  H.  Hsieh. 

2.  "Crack  Propagation  in  Burning  Solid  Propellants,"  presented  in  ONR  Work¬ 
shop  on  Dynamic  Deformation,  Fracture  and  Transient  Combustion,  held  at 
Chestertown,  M),  July  28-29,  1583  by  K.  K.  Kuo. 
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II.  INVESTIGATION  OF  COMPOSITE  PROPELLANT  IGNITION  MODEL 
2 . I  Brief  Introduction  of  Nitramine  Propellant  Ignition 

A  fundamental  understanding  of  the  ignition  process  of  nitramine  solid 
propellants  is  important  to  the  design  of  rocket  motors  and  artillery  propul¬ 
sion  systems.  Most  of  the  solid  propellant  ignition  studies  were  performed 
on  AP-based  composite  propellants  and  homogeneous  propellants.  A  detailed 

review  of  literature  on  solid  propellant  ignition  was  conducted  by  Price  et 

1  2  •  3 

al.,  and  more  recently  by  Hermance  ,  and  Kulkarni,  Kumar,  and  Kuo.  Igni¬ 
tion  of  traditional  AP-based  composite  solid  propellants  located  at  the  tip 
of  an  inert  crack  was  investigated  both  experimentally  and  theoretically  by 
Kuo  and  coworkers. A  comprehensive  model^  and  numerical  solutions  for 
ignition  of  AP-based  composite  solid  propellants  under  rapid  pressurization 
rate  were  presented.  The  rapid  pressurization  situation  is  typical  of  igni¬ 
tion  transients  in  rocket  motors  and  the  onset  of  deflagration-to-detonation 
transition  (DDT)  in  confined,  damaged  propellant  grains. 

Itie  use  of  nitramine  oxidizer  in  composite  propellants  has  become  of 
great  interest  because  it  offers  many  advantages  in  rocket  and  gun  propul¬ 
sion.  High  specific  impulse  is  obtained  due  to  the  combination  of  high 

8  9 

energy  and  low  molecular  weight  of  product  gases.  '  The  excellent  thermal 
stability  of  HMX  and  RDX  makes  nitramine  composite  propellants  the  best  can¬ 
didate  for  Low  Vulnerability  Anmunition  (LOVA)  propellants.^'^  In 
addition,  replacement  of  the  currently  used  AP  oxidizer  by  HMX  or  RDX  oxidiz¬ 
er  eliminates  HC1  in  the  exhaust  gases  and  reduces  secondary  smoke  due  to 
nucleation  of  moisture  to  form  droplets. 


In  the  ignition  of  nitramine  propellants,  experimental  work  conducted  by 

12  13 

DeLuca,  Caveny,  Ohlemiller,  and  Surtmerfield  '  showed  that  HMX  composite 

propellants  are  more  difficult  to  ignite  by  arc-image  furnace  or  CC>2  laser 

than  AP-based  composite  propellants  and  noncatalyzed  double-base  propellants. 

14 

Simile  results  obtained  by  Birk  and  Caveny  also  showed  that  HMX  composite 
propellants  are  more  difficult  to  ignite  under  transient  flew  conditions  than 
single-base,  double-base,  and  triple-base  propellants.  In  spite  of  these 
studies,  the  ignition  process  of  nitramine  propellants  under  rapid  pressuri¬ 
zation  situation  has  not  been  studied  thoroughly. 

2.2  Theoretical  Approach 

2.2.1  Description  of  Physical  Model 

The  physical  processes  that  lead  to  ignition  are  as  follows.  Hot 
combustion  product  gases  from  the  igniter  or  main  chamber  flow  into  the  crack 
cavity  and  pressurize  it.  This  causes  the  pressure  of  the  gases  adjacent  to 
the  propellant  surface  to  increase  rapidly.  As  the  process  continues,  energy 
is  transferred  from  the  hot  gases  to  the  propellant.  Following  a  period  of 
inert  heating,  the  propellant  starts  to  decompose.  The  fuel  and  oxidizer 
species  evolved  from  the  surface  diffuse  and  mix  with  the  surrounding  gases. 
In  the  gas  phase,  oxidizer  species  react  with  the  fuel  species;  at  the  same 
time,  surface  and/or  subsurface  reactions  in  the  solid  phase  continue.  When 
the  net  heat  evolved  from  chemical  reactions  overcomes  the  heat  losses,  the 
temperatures  st£.tt  to  rise.  Eventually,  ignition  characterized  by  a  high 
rate  of  chemical  reactions  and  heat  release,  is  achieved. 

The  physical  model  considers  an  oxidizer  particle  embedded  in  a  fuel 
binder  matrix.  In  an  actual  composite  propellant,  oxidizer  particles  of 
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random  size  (within  a  certain  range)  are  distributed  in  a  fuel  binder. 
Figure  1  shows  typical  side  cross-section  and  top  views  of  a  typical  unimodel 
composite  propellant.  Following  Ref.  7  in  the  mathematical  simulation  of 
this  type  of  propellant,  the  shape  of  the  oxidizer  particles  is  approximated 
by  a  cylindrical  pellet  of  radius  R^,  and  thickness  L,  enveloped  by  a  fuel 
binder  of  radius  R2  as  shown  in  Fig.  2.  A  two-dimensional  geometry  closely 
represents  a  composite  propellant,  vrtiile  being  mathematically  tractable.  The 
gas-phase  pressure-time  (P-t)  relationship  is  assumed  to  be  known  either  from 
a  measured  P-t  trace,  or  a  prescribed  P-t  information  simulating  a  given  op¬ 
erating  condition.  Because  pressure  is  not  usually  extremely  high  during  the 
ignition  interval,  the  perfect  gas  law  can  be  used  as  the  equation  of  state. 


2.2.2  Governing  Equations 

For  the  coordinates  shown  in  Fig.  2,  heat  equations  for  the  oxidizer  and 
fuel  binder  in  the  solid  phase  are 

Oxidizer:  p.  c  H  -  T"  0W  4^)  +  “  ar  I?) 

Ox, a  Ox,s  3t  3z  Ox, a  3z7  r  3r  Ox,s  3r' 


+ 


a"' 

^Ox.s 


»-**  I  ■  A  «■*..  H>  +  T  & *».,«  > 


+  «F„ 


where 


ittl  m  ifll  ^  "III 

HOx,s  s0x,  radiation  q0x, pyrolysis 


(1) 


(2) 


(3) 


1 
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CROSS-SECTIONAL  VIEW 
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OXIDIZER 


FUEL 


an 


A  TYPICAL  ELEMENT 


Fig.  1  Heterogeneous  Structure  o*  Nitramine 
Composite  Propellants 


IGNITER  GAS 


r  (J) 


ADIABATIC 

BOUNDARIES 


Fig.  2  Statistically  Averaged  Element  and  Finite 
Difference  Grid  Pattern 
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pOxz 

^Ox, radiation  *  T0xGS0xe  “  SOxIs 


qOx, pyrolysis  "  "  Z0x,py^0x,py 


Ox,py  u 


The  functional  form  of  pyrolysis  heat  absorption  terra  in  Eq.  (5)  is  the 

same  as  that  used  by  Ohlemiller  and  Sumner  field.^  The  expression  for  q£'s 

is  similar  to  that  of  aJI'  . 

vX  f  s 

The  gas-phase  conservation  equations  are 


^  a 

Continuity  Eq:  -r-®  + 

dt 


'Energy  Eq:  c  p  —  +  p  c  v  ~  ~  »  irl  (k 

P  g  3t  g  P  *  3z  3t  3zi  ^  g  3*' 


+  —  (V  r  + 

+  r  3r  V  3r'  +  qc 


3Yi  3yh  a  !  3*4 

Species  Eq:  Pg  ^  +  p^  -jJ-  -  £  (pg  D \ ^±) 


+  ?3?  V'lT*  (8) 

where  j  *  1,  2,  3,  or  4  stands  for  the  gas-phase  species: 

oxidizer-rich  gases  (called  Oxidizer) ,  N02,  10,  and  fuel-rich  gases  (called 
Fuel) .  The  symbol  10  stands  for  intermediate  oxidizer  to  be  defined  in  a 

»  IV  V  sM  V 

later  section.  The  inhomogeneous  terms  q  and  «  are  obtained  by  using 
chemical  kinetics  information  described  in  a  later  section. 


Since  the  region  of  interest  in  the  gas  phase  during  the  ignition  process 
is  very  small  (in  the  order  of  1  mm),  local  pressure  in  this  region  is  con¬ 
sidered  to  be  spatially  uniform.  However,  the  pressure  is  allowed  to  vary 
with  respect  to  time.  Therefore,  the  gas-phase  momentun  equations  are  re¬ 
placed  by  the  measured 
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P  =  P(t) 


The  equation  of  state  for  the  gas  phase  is 


„  PW 

Mg  R  T 
u 


2.2.3  Initial  and  Boundary  Conditions 


The  initial  condition  for  the  solid  phase  (both  fuel  and  oxidizer)  is 


at  t  *  0  :  ?(0,r,z)  *  (11) 

The  boundary  conditions  for  the  solid  phase  are  as  follows.  The 
temperature  far  from  the  interface  is  the  same  as  the  propellant  initial  tem¬ 
perature,  i.e. 


at  z  -*•  -»  :  T  » 


Temperature  continuity  at  various  interfaces  requires 


r  -  R 


T  .  T 

+  "  _  _  „  -  and 

l  r  "  Kj_  z  -  -L 


-  -L+ 


On  symmetric  surfaces,  adiabatic  conditions  are 


on  r  *  0 


and  on  r  >  R. 


At  solid-gas  and  oxidizer-fuel  interfaces,  heat  flux  balance  gives 

ats-v8  *  fS| .  „  ♦  v,«  ♦ 
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r.  p  T(c  -  c  )  +  q" 
bs  8  ps  pg  s-g 

where  q^  is  the  net  heat  generation  at  the  surface. 


(16) 


at  z  «  -L 


3T 


•  k  — 

*  *7t«  3z 


-k  il 

-L  Ox, a  3z 


-L+  + 


-a  l 

Ox 


a0x-FT0xG8 


at  r  ■  R^  and  -L  <  z  <  0 


(17) 


.  .  t  3T 

* »*  3r  Rj  + 


-  k.  & 


(18) 


Ox,*  3r|  R^* 

Initial  conditions  for  the  gas-phase  equations  (for  j*l,2,3,4)  are: 

at  t  »  0  :  v,(0,r,z)  ■  0  ;  T(0,r,z)  -  T  Y,(0,r,z)  -  Y.  . 

*■  g,i  j  J*4, 

(19) 

Boundary  conditions  for  the  gas  phase  are  as  follows.  The  syrmetry 
conditions  at  the  centerline  and  outer  boundary  give 
3T 


3Y 


on  r  ■  0 


3r 


-  0  ; 


3r 


(20) 


3T  3Yi 

and  on  r  ■  R2  :  37  “  0  »  ■  0 


(21) 


Far  away  from  the  surface 
3T  3Yi 

z  -*■  o«  ;  2JL  m  0  •  — J-  -  0 

3z  3z 


(22) 


The  temperature  continuity  at  the  interface  gives 


on  z  ■  z 


a-g 


-  T  | 


-s-g+  ' z*-g- 


g 


(23) 


The  overall  and  individual  species  mass-flux  balance  at  the  solid-gas  inter¬ 
face  is  given  respectively  by 


on  z  «  z 


s-g 


PgVz  "  psrb 


(24) 


-  / 


■/ 
/  • 


\  / 


•\l  <•.  ft  A  V.V.V: V.V  V  V  ">  *„-  *.»  *>  -  •  -  .•>  -  •  ^  >  •/.  -  V  ' 
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*  p 

This  set  of  partial  differential  equations  is  being  solved  numerically  and 
the  solution  will  be  presented  in  the  future. 

2.2.4  Chemical  Kinetics  of  Nitramine  Composite  Propellants 

It  is  generally  known  that  the  decomposition  of  nitramine  dominates  the 
overall  chemical  process  in  deflagration  of  nitramine  composite  propellants. 
A  proper  starting  point  in  understanding  the  chemical  kinetics  of  combustion 
of  nitramine  composite  propellants  is  to  investigate  the  pure  nitramine  de¬ 
composition.  Due  to  their  complexity  and  special  characteristics,  the 
decomposition  kinetics  of  nitramines  have  been  investigated  extensively 

during  the  last  thirty  years.  Recently,  Schroeder1®-1®,  Boggs and 
20 

McCarty  have  published  comprehensive  review  articles  on  the  kinetics  of  ni- 
tr amine.  Schroeder  and  Boggs  studied  the  subject  from  an  organic  chemist's 
point  of  view  while  McCarty  approached  it  from  a  combustion  engineer's  point 
of  view.  Primary  N-NC>2  cleavage  is  regarded  as  the  important  initial  decom- 
position  step,  but  there  appears  to  be  some  evidence  supporting  two  other 
possible  initial  steps  in  the  thermal  decomposition  of  HMX  and  RDX:  these  are 

unimolecular  elimination  reaction  giving  HONO  and  C=N  double  bond  in  the  ring 

22  23 

instead  of  CH^jNOj'  and  unimolecular  C-N  bond  cleavage.  Brill  and  Kar- 

24 

powicz  suggested  that  the  initial  rate  controlling  step  for  decomposition 
in  the  condensed  phase  involves  mainly  the  disruption  of  the  strong  intermo- 
lecular  electro-static  forces  between  HMX  molecules,  and  also  between  the  HMX 
molecules  and  its  decomposition  products.  Their  postulation  is  based  upon  a 
comparison  of  Arrhenius  data  for  phase  transformation  with  data  reported  for 
the  condensed  phase  decomposition  process.  The  activation  energies  of  the 


and  p  v  Y. 

&  z  j 


v  Y . 
8  *  j 


Z  4. 
S-g+ 


3Y 

V  -r-i  - 
3z 
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phase  transitions  reflect  the  energy  needed  to  disrupt  electro-static  forces 
in  the  crystal  lattice  of  HMX. 


r  : 


Follow-up  reaction  steps  can  be  divided  into  two  categories:  unimolecular 

23 

and  bimolecular .  It  appears  that  stepwise  loss  of  CH2NNO2  rosy  he  an  impor¬ 
tant  unimolecular  follow-up  reaction,  and  that  bimolecular  reactions  of 

radical  or  products  formed  in  the  decomposition  with  unreacted  HMX  or  RDX  mo- 

23 

lecules  is  also  important.  Suryanarayana  et  al.  proposed  a  concerted 

intramolecular  condensed  phase  decomposition  mechanism  initiated  by  C-N  bond 

25  26 

fission.  Gosgrove  and  Owen  '  suggested  a  gas  phase  mechanism  initiated  by 

27 

N-N  bond  fission.  Goshganan  measured  detailed  concentration  and  tempera¬ 
ture  profiles  for  RDX  and  HMX  by  using  mass  spectrometry.  The  presence  of 
NOj  under  both  molten  and  deflagration  condition  was  indicated,  verifying  the 
N-N  cleavage  is  the  important  initial  step  for  decomposition,  and  the  N02  is 
a  reactant  of  the  secondary  reactions  in  the  vapor  phase.  Unfortunately, 
most  investigations  of  deconposition  kinetics  of  nit  rami  nes  were  conducted  at 
low-tenperature  and  low-pressure  situations  which  may  not  simulate  deflagr¬ 
ation  of  these  substances  in  rocket  motors.  Detailed  mechanisms  obtained  by 
organic  chemists  are  too  complex  to  apply  to  actual  highly  transient  ignition 
situations. 


For  kinetic  mechanism  in  a  combustion  model,  Ben-Reuven  and  Caveny28  pro¬ 


posed  a  comprehensive  kinetic  model  for  HMX  as  the  monopropellant,  by 
considering  a  reacting  molten  layer  coupled  to  the  gas  field  through  conserve- 
ation  conditions  satisfied  by  all  chemical  species  and  enthalpy.  The  primary 
initial  overall  reaction  used  in  Ben-Reuven  and  Caveny 1  s  gas-phase  model  of 
deflagration  of  HMX  is  assumed  to  be  similar  to  RDX,  and  is  represented  by 


-* -.'.'VY'- -*  V.V.7  .V;.',v.r5r;  Y*.nr-^r'  J-  yjy  I 
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HMX(g) - -  4/3N02  +  4/3N20  +  2N2  +  4CH2<)  (26) 

Since  nitrogen  dioxide  reacts  quite  rapidly  with  formaldehyde,  the  overall 
consecutive  reaction  in  the  gas  phase  is  represented  by 

7N02  +  5CH20  - -  7N0  +  3C0  +  2C02  +  SHjO  (27) 

In  the  liquid  phase,  the  overall  reaction  step  was  represented  by 
HMX(1) - -  4CH20  +  3N20  +  NO  +  1/2N2  (28) 


For  RDX,  a  similar  kinetic  scheme  was  used  by  Ben-Reuven  and  Summer - 


field29  and 

is  represented  by: 

i 

t 

““(!)  — 

♦ 3CH^0  +  9/4N20  +  3/4NO  +  3/8N£ 

(29) 

““(g)  — 

-  N02  |+  N20  f  3CH20  +  3/2N2 

(30) 

the  secondary  overall  reaction  is  also  the  oxidization  of  CH_0  by  NO.  and  is 

i  c 

\ 

not  repeated  here1. 


In  the  ignition-  and  combustion  study  of  nitramine  composite  propellants, 

i  ■ 

most  attention  in  kinetic  mechanism  has  been  devoted  to  nitramine  fillers; 
very  little  attention  has  been  given  to  the  binder.  However,  in  order  to 
gain  deeper  understanding  of  the  ignition  process  and  to  improve  the  pred¬ 
ictability  of  theoretical  model,  it  is  necessary  to  consider  the  role  of  the 

binder  in  the  ignition  and  combustion  process.  Experimental  results  of 
30 

Kubota  et  al.  revealed  that  flame  structure,  combustion  mechanism,  and 
burning  rate  of  the  HMX  composite  propellant  are  strongly  influenced  by  the 
type  of  binders  used.  When  a  high-oxygen  content  binder  (e.g.. 


HTPS,  HTPE, 
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or  HTPA)  was  used,  flame  structures  became  homogeneous  because  of  the  form¬ 
ation  of  a  homogeneously  mixed  HMX/binder  molten  layer  on  the  burning 
surface.  On  the  other  hand,  when  a  low  oxygen  content  binder  such  as  HTPB 
was  used,  a  large  amount  of  carbonaceous  material  was  formed  on  the  burning 
surface,  making  the  structures  of  the  burning  surface  and  the  gas  phase 
highly  heterogeneous.  The  decomposition  reactions  of  polymeric  binder  are 
extremely  complicated  and  are  not  yet  fully  understood.  Hie  primary  decom¬ 
position  products  are  possibly  aldehydes,  hydrocarbons,  and  hydrogen.  The 
gaseous  mixture  of  the  decomposed  gases  of  the  binder  is  ejected  into  the  gas 
phase  from  the  burning  surface  of  propellant  and  mixed  again  with  gases  de¬ 
composed  from  nitramine  oxidizers.  Hie  mixture  exhibits  rapid  exothermic 
reactions  which  include  many  parallel  and  successive  reactions. 

In  modeling  the  kinetics  of  nitramine  composite  propellant  under  highly 
transient  ignition  conditions,  the  basic  idea  follows  the  overall  reaction 
representation  of  Ben-Reir'en  and  Caveny  in  order  to  avoid  complexity  of  de¬ 
tailed  kinetic  mechanism.  However,  heterogeneous  structure  rather  than 
homogeneous  molten  layer  is  still  considered  valid,  before  onset  of  steady 
deflagration  of  nitramine  composite  propellants.  Hiis  is  consistent  with  the 
experimental  test  conditions  which  are  highly  transient  with  the  enormous 
heating  rate  at  the  propellant  surface.  A  chemical  reaction  mechanism  of  the 
ignition  of  nitramine  composite  propellant  is  proposed.  It  consists  of  the 
following  reactions: 

(a)  following  the  vaporization  of  nitramine  crystals,  primary  gas-phase 
reaction  occurs  due  to  decomposition  of  HMX  (or  RDX)  to  form  NC>2  and 
Intermediate  Oxidizers  (10)  which  consists  of  other  species  than  N02 


■■j  r.  it;  r 


'jMji  ■  ■ .  u  mjwyvjn 


:  wr.  tv  rji 
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in  Eq.  (26)  or  (30) ; 

(b)  reaction  of  10  species  with  NC>2  t0  fonn  oxidizer  gases  (this  reaction 
is  considered  to  be  the  secondary  gas-phase  reaction) ; 

(c)  decomposition  and  pyrolysis  of  binder  into  fuel  gas; 

(d)  oxidation  of  fuel  gases  by  oxidizer  gases  to  form  final  product 
gases;  and 

(e)  heterogeneous  reaction  between  oxidizer  gases  and  solid  fuel  binder 
on  the  propellant  surface  to  form  final  product  gas. 

The  proposed  simplified  kinetic  schemes  and  data  for  three  types  of  ni- 
tramine  composite  propellants  considered  in  this  paper  are  given  in  Tables 
1-3.  The  kinetic  data  for  reaction  (d)  in  these  three  tables  are  based  upon 
the  representative  species  which  have  the  highest  concentrations  in  oxidizer 
and  fuel  gases.  In  most  cases,  fuel  gases  are  represented  by  CH^,  and  oxi¬ 
dizer  gases  are  represented  by  NO  and  NjO.31  Kinetic  data  of  reactions  (e) 
are  largely  unavailable.  Representative  data  may  be  obtained  in  the  future 
for  major  species  involved  in  heterogeneous  reactions. 

2.3  Experimental  Approach 

2.3.1  Test  Apparatus 

The  main  test  chamber  with  the  test  configuration  was  designed  to  be 
pressurized  by  high-temperature  and  high-pressure  gases  generated  in  a  solid- 
propellant  driving  motor.  Pressurization  rates  of  the  test  chamber  closely 
simulate  conditions  during  the  ignition  transient  of  actual  rocket  motors. 
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Table  1.  Kinetic  Scheme  and  Data  for  Ignition  of 
HMX/CTPB  Composite  Propellants 


^^^-^Kinetic  Data 

Reaction 

Pre¬ 

exponential 

Factor 

Heat  of 
Reaction 

Source 

(a)  HMX-»-N02  +  7.67  10 

IxlO16'4 

1/s 

1.931x10® 

J/kaole 

-5.166x10® 

J/kaole 

Ben-Reuven 

Caveny 

(1979) 

(b)  N02  +  7.67  10 -►  9.62  OX 

lxlO6 

m2/kmole-s 

0.802x10® 

J/kmole 

-1.923x10® 

J/kmole 

Ben-Reuven 

Caveny 

(1979) 

(c)  CTPB  -*■  F 

128 

kg/m2-s 

0.4389x10® 

J/kmole 

1.59xl06 

J/kg 

Cohen 

Fleming 

Derr 

(1974) 

(d)  4.361  OX  +  F-^8.668  P 

lxlO16 

mJ/kmole-s 

0.4x10® 

J/kmole 

-7.4xl07 

J/kmole 

Drummond* 

(1969) 

(e)  OX  +  CTPB-*-1.99  P 

— 

~ - 

— 

Kinetic  data  are  based  on  the  representative  species  which  have  the  highest 
concentration  in  oxidizer  and  fuel  gases. 


/ 
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Table  2.  Kinetic  Scheme  and  Data  for  Ignition  of 
RDX/CTPB  Composite  Propellants 


^^x^Kinetic  Data 

Reaction 

Pre¬ 

exponential 

Factor 

Activation 

Energy 

Heat  of 
Reaction 

Source 

(a)  RDX-*-N02  +  5.5  10 

lxlO15,5 

1/s 

1. 73x10® 
J/kmole 

-3.97x10® 

J/kmole 

Ben-Reuven 

Caveny 

Summer field 

Vichnevetsky 

(1976) 

lxlO9 

m^/kmole-s 

0.79x10® 

J/kmole 

-1. 923xl08 
J/kmole 

Ben-Reuven 

Caveny 

Vichnevetsky 

Sunmerfleld 

(1976) 

(c)  CTPB— ►  F 

128 

kg/m2-s 

• 

0.4389x10® 

J/kmole 

1.59xl06 

J/kg 

Cohen . 

Fleming 

Derr 

(1974) 

(d)  3.71  OX  +  F->- 8.668  P 

0.95xl016 

nr/kmole-s 

0.4xl08 

J/kmole 

-7.2xl07 

J/kmole 

Drummond* 

(1969) 

1 

H 

i 

i 

i 

i 

1 

i 

i 

i 

1 

— 

Kinetic  data  are  based  on  the  representative  species  which  have  the 
highest  concentration  in  oxidizer  and  fuel  gases. 
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Table  3.  Kinetic  Scheme  and  Data  for  Ignition  of 
HMX/HTPB  Composite  Propellants 


^^^^'^^jCinetic  Data 
Reaction 

Pre¬ 

exponential 

Factor 

Activation 

Energy 

Heat  of 
Reaction 

Source 

(a)  HMX-*‘N02  +  7.67  10 

IxlO16-4 

1/s 

1. 931xl08 
J/kmole 

1 

-5.166xl08 

J/kmole 

Ben-Reuven 

Caveny 

(1979) 

(b)  N02  +  7.67  10-*-  9.62  OX 

lxlO6 

m^/kmole-s 

0.802xl08 

J/kmole 

8 

-1.923x10° 

J/kmole 

Ben-Reuven 

Caveny 

(1979) 

(c)  HTPB  -*F 

2990 

kg/m2-s 

0. 7064xl08 
J/kmole 

1.81xl06 

J/kg 

Cohen 

Fleming 

Derr 

(1974) 

(d)  4.290  OX  +  F  -*  8.53  P 

0. 9xl016 
nr/kmole-s 

0. 39xl08 
J/kmole 

-7.3xl07 

J/kmole 

Drummond* 

(1969) 

(e)  OX  +  HTPB  -*1.99  P 

— 

— 

— 

— 

Kinetic  data  are  based  on  t’ e  representative  species  which  have  the  highest 
concentration  in  oxidizer  and  fuel  gases. 
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High  pressurization  rates  obtained  in  this  setup  can  also  simulate  those  con¬ 
ditions  encountered  during  onset  of  DDT.  Using  this  test  rig,  it  is  possible 
to  obtain  pressures  >ip  to  50  MPa  ( ~  70CX3  psi),  and  pressurization  rates  up  to 
150  GPa/s  ('1.5  x  106  atm/s). 

A  schematic  diagram  of  the  test  section  is  shown  in  Fig.  3.  A  crack-like 
cavity,  formed  between  a  transite  slab  and  a  plexiglass  window,  is  situated 
normal  to  the  flow  direction  of  gases  in  the  main  chamber.  The  gap  width  of 
the  crack  is  determined  by  the  recessed  depth  of  the  transite  slab  in  the 
test  configuration.  Because  the  main  motivation  of  this  work  is  to  study  the 
effect  of  propellant  ingredients,  the  gap  width  was  held  constant,  and  equal 
to  1.27  ires.  The  cavity  length  for  the  test  configuration  was  161  mm.  As 
shown  in  Fig.  3,  two  propellant  test  samples  are  located  at  the  tip  of  the 
cavity.  Two  pressure  transducer  ports  are  provided  in  the  crack  cavity:  one 
near  the  test  sample  shown  in  Fig.  3,  and  the  other  at  the  crack  entrance. 
The  ignition  event  is  observed  through  a  set  of  transparent  plexiglass  win¬ 
dows.  The  window  assembly  is  held  in  place  by  a  window  retainer,  which  forms 
the  front  half  of  the  chamber  shown  in  Fig.  4.  With  the  exception  of  the 
exit  nozzle  through  which  gases  are  discharged  into  the  atmosphere,  the  cham¬ 
ber  is  completely  sealed  during  tests.  For  safe  operation,  the  chamber  is 
also  equipped  with  a  safety  head,  which  houses  a  bursting  diaphragm;  the  dia¬ 
phragm  ruptures  at  a  prespecified  pressure  of  70  MPa  (-  10,000  psi).  The 
driving  motor  system  used  in  this  study  is  a  small  solid-propellant  motor  in¬ 
itiated  by  an  electric  primer.  Details  of  the  driving  motor  system  can  be 
found  in  Ref.  32.  The  high-temperature  combustion  product  gases  generated  in 
this  motor  flow  through  a  multiperf orated  nozzle  into  the  main  test  chamber, 
and  pressurize  the  chamber  at  high  loading  rate. 
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Fig.  4  Photograph  of  the  Assembled  Test  Rig 


It 


22 


Figure  4  shows  the  assembled  test  chamber.  The  driving  motor  is  located 
inmediately  below  the  base  mounting  plate,  and  hence  cannot  be  seen  in  this 
photograph.  Exhaust  gases  from  the  driving  motor  flow  vertically  upward  into 
the  main  flow  channel  located  at  the  left  end  of  the  test  chamber.  The  test 
chamber  has  two  window  openings.  In  order  to  have  a  close-up  view  of  the 
ignition  process  near  the  right  most  portion  of  the  longer  window,  a  light 
pipe  and  a  prism  were  used  to  give  light  signals  at  the  onset  of  the  hot  gas 
flow  in  the  test  diamber.  By  so  doing,  the  event  of  onset  of  hot  gases  (sig¬ 
nified  by  a  bright  circular  spot)  and  the  ignition  process  can  be 
photographed  simultaneously  on  the  high-speed  film.  Before  most  test  fir¬ 
ings,  a  nitrogen  feedline  is  turned  on  to  purge  the  chamber,  eliminating 
oxygen  in  the  test  chamber  during  initial  ignition.  Purged  gases  leave  the 
test  chamber  through  the  exhaust  pipe  mounted  vertically  on  the  left  side  of 
the  test  chamber.  The  purpose  of  driving  oxygen  from  the  system  is  to  reduce 
the  complexity  of  the  ignition  process  so  that  oxygen  will  not  participate  in 
the  ignition  and  combustion  processes. 

Special  features  of  the  test  configuration  are  near-infrared  photodiodes 
used  to  sense  infrared  radiation  from  surfaces  of  propellant  samples,  and  the 
sample  holder  configuration  is  designed  to  provide  a  clear  view  of  the  igni¬ 
tion  event.  The  near-infrared  photodiode  is  the  product  of  Texas  Instruments 
Company  (T.  I.  No.  1N2175) ,  with  photocurrent  rise  time  about  2  ys  and  spec¬ 
tral  response  range  from  0.42  ym  to  1.12  ym.  As  shown  in  Fig.  5,  two 
near-infrared  photodiodes  are  installed  into  the  chamber.  Infrared  light  em¬ 
itted  from  the  burning  surface  of  the  propellant  sample  passes  through  a 
fiber  optic  rod  and  a  glass  rod  window,  and  is  sensed  by  the  near-infrared 
photodiode.  To  measure  the  sensitivity  of  the  near-infrared  photodiode  in 
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Fig.  5  Detailed  Schematic  Diagram  of 

the  Setup  in  Test  Configuration  B 
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response  to  gas-phase  radiation  from  propellant  flames,  separate  tests  were 
conducted.  It  was  found  that  the  gas-phase  radiation  cannot  be  detected  by 
the  near-infrared  photodiode;  only  the  near-infrared  light  from  the  burning 
solid  surface  can  be  monitored.  This  implies  that  the  signal  recorded  by  the 
sensor  represents  the  surface  temperature  rise  rather  than  the  arrival  of  hot 
gases.  Based  upon  the  distinct  feature  of  the  near-infrared  sensor,  this 
test  configuration  was  used  in  most  of  the  test  firings.  An  external  photo¬ 
diode  system,  consisting  of  photographic  lens  mounted  on  an  optical  bed,  a 
fast-response  photodiode,  and  an  amplifier,  was  also  used  in  the  test  confi¬ 
guration. 


The  data  acquisition  system  used  in  the  present  study  is  very  similar  to 
that  described  in  the  paper  by  Siefert  and  Kuo.  This  system  is  comprised 
of  (1)  pressure  measuring  system,  (2)  near-infrared  photodiode  (or  ionization 
pin)  setup,  (3)  high-speed  photography  system,  (4)  photodiode  system,  and  (5) 
transient  waveform  recording  system.  Piezoelectric  quartz  transducers  with  a 
rise  time  of  1.5  ys  and  natural  frequency  of  300  kHz  were  used  to  measure  the 
pressure.  The  near-infrared  photodiode  system,  photodiode  system,  and  a  16 
mm  Hycam  movie  camera  were  used  in  most  of  the  tests,  and  the  time  interval 
between  consecutive  pictures  is  about  60  ys.  At  ignition,  it  is  expected 
that  j  the  near-infrared  photodiode  signal  will  begin  to  rise  abruptly,  and 
that  a  bright  light  will  appear  at  the  propellant  surface  on  the  film. 
Therefore,  to  determine  the  onset  of  ignition,  the  following  procedure,  in 
which  the  film  is  read  in  conjunction  with  the  photodiode  system,  is  adopted: 

(1)  the  picture  that  shows  the  first  occurrence  of  bright  light  is  de- 

’ 

fined  as  tfiM, 


25 


(2)  the  corresponding  point  is  located  on  the  photodiode  signal 

(3)  the  point  of  abrupt  rise  in  the  photodiode  signal  within 
tfi^  is  taxen  as  the  actual  time  for  onset  of  ignition. 


2.3.2  Propellants  Used  in  Experiments 


,  and 

30  Ps  of 


Three  types  of  nitramine  propellants  (1,  2,  and  3)  were  used  in  experi¬ 
mental  study.  Propellant  1  is  HMX-based  composite  propellant  with  79%  HMX 
and  21%  CTPB  binder  mix.  Propellant  2  is  RDX-based  composite  propellant  with 
79%  RDX  and  21%  CTPB  binder  mix.  The  comparison  of  propellant  1  and  2  gives 
the  effect  of  nitramine  filler.  Propellant  3  is  HMX-based  composite  propel¬ 
lant  with  80%  HMX  and  20%  HTPB  binder  mix.  The  binder  effect  can  also  be 
studied  by  comparing  the  results  of  propellant  1  and  3.  Several  important 
thermal  properties  of  propellants  and  their  products  are  tabulated  in  Table 
4.  Some  of  the  physical  properties  of  HMX,  RDX,  HTPB,  and  CTPB  are  tabulated 
in  Table  5. 

2.4  Results  and  Discussion 

Figure  6  shows  a  typical  set  of  time  correlated  data  which  includes 
trigger  pulse  (common  time  signal) ,  pressure  at  the  crack  entrance  (G2) , 
pressure  at  the  crack  tip  (G3) ,  near-infrared  photodiode  signals  (Nil  and 
NI2)  for  propellant  samples  HMX/CTPB  and  HMX/HTPB,  amid  photodiode  signal  in 
the  region  near  HMX/CTPB  sample.  Trigger  pulse  is  generated  by  a  light- 
emitting  diode  (LED)  driver  unit,  and  is  used  to  correlate  time  between  data 
obtained  from  high-speed  movie  films  and  transient  waveform  recording 
devices.  When  hot  gases  from  the  driving  motor  reach  the  crack  entrance, 
pressure  at  the  entrance  increases  rapidly.  The  first  discernible  pressure 
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Table  4.  Thermal  Properties  of  Three  Types  of  Nitramine 
Composite  Propellants  and  Their  Products 


Av.  Mol.  Wt.  of  Product  Gases 
(kg/kmole) 


19.59 


19.58 


19.33 


Table  5.  Properties  of  Nitramine  Propellant  Fillers  and  Binders 
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rise  at  the  tip  occurs  some  time  later  than  the  corresponding  pressure  rise 
at  the  crack  entrance.  The  rising  portion  of  the  pressure  trace  at  the  crack 
tip  is  quite  linear,  and  the  averaged  pressurization  rate  at  the  tip  can  be 
obtained  easily.  In  order  to  measure  the  time  required  for  the  onset  of  gas 
evolution,  the  first  discernible  pressure  rise  at  the  tip  is  chosen  as  the 
time  equal  to  zero.  Shortly  after  the  pressure  rise  at  the  tip,  the  near- 
infrared  photodiode  signals  for  the  propellant  samples  show  abrupt  rises,  es¬ 
pecially  for  Nil.  These  rises  correspond  to  the  increase  of  the  surface 
temperatures  of  propellant  samples.  After  analyzing  the  motion  picture 
films,-  the  onsets  of  gas  evolution  of  two  propellant  samples  are  marked  as 
two  black  dots  on  the  near-infrared  photodiode  signals.  As  shown  in  Fig.  6, 
the  onset  of  gas  evolution  occurs  during  the  uprising  portion  of  the  near- 
infrared  photodiode  signal,  rather  than  the  peak  of  the  near-infrared  photo¬ 
diode  signal.  The  external  photodiode  signal  (PD)  for  detecting  light 
intensity  near  the  HMX/CTPB  propellant  sample  also  shows  a  similar  rise  with 
Nil. 

The  consistency  between  the  two  independent  measuring  systems  is  helpful 
in  determining  the  time  of  onset  of  gas  evolution.  In  this  particular  case, 
the  HMX/CTPB  propellant  begins  to  gasify  slightly  sooner  than  the  HMX/HTPB 
propellant.  After  momentary  ignition,  both  propellant  samples  were  extin¬ 
guished;  thus  the  two  black  dots  in  Fig.  6  can  be  regarded  only  as  onset  of 
gas  evolution,  and  not  as  onset  of  ignition.  After  the  strong  pulse,  near- 
infrared  photodiode  signals  decrease,  due  to  the  extinguishing  of  propellant 
samples,  and  in  part  to  the  coating  on  the  end  of  the  fiber  optic  by  product 
gases  fran  the  driving  motor  and  pyrolyzed  gases  from  propellant  samples. 


G2  (Crack  Entrance), 


0123  100  200  300 

Time,  t,  ms 

Fig.  6  A  Typical  Set  of  Time  Correlated  Data  of  a 
Test  Firing  with  Two  Different  HMX-Based 
Composite  Propellants 
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Figure  7  shows  the  experimental  data  of  onset  of  gas  evolution  of  two 
HMX-based  propellant  samples  as  a  function  of  the  pressurization  rate.  These 
two  propellants  have  nearly  the  same  weight  percentage  of  HMX,  but  with  dif¬ 
ferent  binder  mixes  (CTPB  and  HTPB)  .  Although  a  few  tests  using  these  -..wo 
propellants  under  identical  operating  conditions  suggest  that  the  time  re¬ 
quired  for  onset  of  gas  evolution  is  slightly  shorter  for  HMX/CTPB 
propellant,  there  is  no  absolute  indication  that  this  is  true  for  all  circum¬ 
stances. 

Data  shown  in  Fig.  7  were  compiled  from  a  number  of  test  firings,  and 
could  be  considered  more  or  less  one  family  and  fitted  by  a  single  curve. 
One  possible  reason  for  this  is  that  the  low  percentage  of  weight  fraction  of 
binder  (  -  201)  does  not  make  the  binder  effect  very  pronounced.  It  is  quite 
obvious  that  the  time  required  for  onset  of  gas  evolution  decreases  with  the 
increase  of  pressurization  rate.  It  should  be  noted  that  none  of  these  data 
corresponds  to  sustained  ignition;  hence  the  curve  can  only  be  regarded  as 
the  boundary  of  onset  of  gas  evolution.  One  major  reason  for  the  fact  that 
the  tests  did  not  sustain  ignition  is  the  short  action  time  of  the  solid  pro¬ 
pellant  driving  motor  (less  than  50  ms). 

The  weak  effect  of  binder  may  be  explained  by  the  following 
consideration.  As  it  is  noted  above  that  in  several  test  firings  conducted 
with  HMX/CTPB  and  HMX/HTFB  as  samples  in  the  same  test  firing,  under  the  same 
operating  situation,  the  time  required  for  onset  of  gas  evolution  of  HMX/HTPB 
was  slightly  longer  than  that  of  HMX/CTPB.  Table  5  shows  that  the  thermal 
diffusivity  of  CTPB  is  larger  than  that  of  HTPB,  which  implies  that  under  the 
same  heating  condition  the  surface  temperature  of  HTPB  will  rise  faster  than 


that  of  CTPB.  However,  using  Table  1  and  3,  it  can  be  seen  that  the  activ¬ 
ation  energy  of  the  pyrolysis  of  CTPB  is  less  than  that  of  HTPB,  which 
implies  that  the  energy  barrier  of  pyrolysis  of  CTPB  is  lower  than  that  of 
HTPB.  Although  the  thermal  diffusivity  effect  and  activation  energy  effect 
work  in  opposite  directions,  the  experimental  data  show  that  the  onset  gas 
evolution  of  HMX/CTPB  is  earlier  than  that  of  HMX/HTPB.  Thus,  it  is  reason¬ 
able  to  conclude  that  activation  energy  of  pyrolysis  of  the  binder  is  more 
important  than  the  thermal  diffusivity  effect. 

Figure  8  shows  the  comparison  of  the  times  required  for  gas  evolution  of 
RDX/CTPB  propellant  (dashed  line)  and  HMX-based  composite  propellants  (shown 
as  solid  line  transferred  from  Fig.  7) .  On  the  same  figure,  the  ignition 
delay  of  AP-based  composite  propellant  is  also  plotted,  in  this  figure,  all 
data  of  AP-based  composite  propellants  correspond  to  sustained  burning  after 
ignition,  while  only  one  datum  of  RDX/CTPB  (marked  by  *)  at  the  lowest  dp/dt 
corresponds  to  sustained  burning.  It  is  well  known,  that  the  difference 
between  the  time  required  for  gas  evolution  and  ignition  is  a  function  of 
heating  rate.  This  difference  is  very  pronounced  for  nitramine  propellants. 
The  time  difference  becomes  smaller  at  low  heating  rates  (or  dp/dt).  Thus, 
it  is  not  surprising  to  obtain  sustained  burning  of  RDX/CTPB  at  low  dp/dt. 
It  is  shown  in  this  comparison  that  the  ease  of  ignitability  increases  in  the 
following  order:  HMX-based,  RDX-based,  and  AP-based  composite  propellants. 
From  Tables  1  and  2,  one  can  see  that  the  activation  energy  of  decomposition 
of  RDX  is  lower  than  that  of  HMX.  Also,  Table  4  shows  that  the  thermal  dif¬ 
fusivity  of  RDX/CTPB  propellant  is  lower  than  that  of  HMX/CTPB.  Therefore, 
RDX/CTPB  is  easier  to  reach  onset  of  gas  evolution  than  HMX/CTPB.  Obvious 
divergences  of  three  boundary  curves  of  HMX-based,  RDX-based,  and  AP-based 
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composite  propellants  can  be  observed  at  low  dp/dt,  while  these  divergences 
tend  to  diminish  at  high  dp/dt.  The  observed  convergence  of  these  boundaries 
at  high  pressurization  rates  is  mostly  due  to  the  high  heating  rates  which 
effectively  reduce  the  influence  of  both  thermal  diffusivity  and  chemical  ki¬ 
netics. 

Figure  9  shows  a  series  of  photographs  developed  from  a  high-speed  movie 
film  obtained  in  a  test  firing  using  the  test  configuration.  This  photograph 
sequence  is  helpful  in  illustrating  the  ignition  event.  The  crack  tip  is  lo¬ 
cated  on  the  right  hand  side  of  the  photograph.  Hot  gases  generated  in 
driving  motor  travel  from  left  to  rights  At  the  lower  left  corner  of  each 

i 

photograph,  there  is  a  prism  which  shows  the  onset  of  light  through  a  light 
pipe  (see  Fig.  4)  upon  the  arrival  of  the  hot  gases  at  the  crack  entrance. 
In  this  test  firing,  RDX/CTPB  propellant  sample  is  installed  at  the  lower 
portion  and  HMX/CTPB  propellant  sample  is  located  at  the  upper  portion.  The 

I 

crack  tip  region  is  pressurized  at|  a  rate  of  16  GPa/s.  The  time  period 
between  consecutive  pictures  is  about  65.6  us.  Picture  No.  1  shows  a  dim  dot 
of  light  on  the  prism  indicating  the  arrival  of  hot  gases  at  the  crack  en¬ 
trance.  After  the  reflection  of  the  compression  wave  at  the  crack  tip,  a 
bright  region  due  to  further  reactions  of  hot  gases  is  initiated  at  the  crack 
tip  (Picture  NO.  2).  The  reaction  front  then  propagates  to  the  left  at  a  ve¬ 
locity  of  88  m/s  relative  to  the  stationary  coordinate.  This  reaction  front 
propagation  causes  the  expansion  of  the  bright  gas  region  in  Pictures  3-8. 
As  time  goes  on,  the  luminous  zone  is  reduced;  this  is  probably  due  to  the 
continuous  chamber  pressurization  by  product  gases  from  the  driving  motor. 
The  incoming  gases  driven  by  the  stronger  pressure  gradient,  compress  the 
gas-phase  reaction  region  and  reduce  the  size  of  the  luminous  zone.  The  re- 
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duction  of  luminous  zone  is  also  partly  caused  by  the  local  pressure  increase 
which  shortens  the  reaction  time.  These  phenomena  are  shown  in  Pictures  8- 
15.  In  Picture  11,  the  surface  of  RDX/CTPB  propellant  sample  starts  to  show 
some  brightness,  which  indicates  the  onset  of  gas  evolution  of  RDX/CTPB  pro¬ 
pellant,  while  the  surface  of  HMX/CTPB  propellant  is  still  dark.  [It  should 
be  noted  that  Pictures  8,  9,  and  10  also  have  some  bright  lines  near  the 
crack  tip;  however,  these  bright  lines  are  due  to  the  reflection  of  light 
from  metal  surfaces  of  the  sample  holder.]  Picture  12  shows  that  HMX/CTPB 
sample  begins  to  gasify.  The  HMX/CTPB  sample  remains  ignited  until  Picture 
23  and  then  extinguishs  thereafter.  The  RDX/CTPB  sample  is  nearly  extin¬ 
guished  as  shown  in  Picture  23  and  it  becomes  definitely  ignited  again  as 
shown  in  Picture  24.  This  ignition  was  sustained  throughout  the  test  firing. 
In  the  later  phase,  more  bright  gases  from  the  driving  motor  enter  the  crack 
cavity.  Due  to  the  increase  in  luminosity,  and  slight  gas  penetration 
between  the  plexiglass  window  and  protruded  surfaces  of  the  transite  base 
plate,  the  screw  heads  on  both  sides  of  the  channel  become  visible  as  black 
circles.  It  is  interesting  to  note  that  the  last  ten  pictures  in  Fig.  9 
showed  some  burning  particles  above  the  RDX/CTPB  surface.  These  particles 
could  either  come  from  the  burning  propellant  sample  due  to  ejection  of 
micro-mechanical  fractured  pieces  or  from  the  driving  motor  system.  At  the 
present  time,  the  source  has  not  been  fully  identified,  and  this  phenomenon 
is  subjected  for  further  investigation. 

The  momentary-ignition  (gas  evolution  processes  followed  by  extinction) 
phenomenon  of  HMX/CTPB  propellant  can  be  explained  as  follows.  When  heating 
rate  is  very  high,  a  thin  thermal  layer  near  the  propellant  surface  is  gener¬ 
ated.  If  the  burning  rate  of  this  hot  thin  layer  of  the  propellant  is  higher 


37 


8 


s' 


£ 

a 

v. 

ft 


<•, 

A 

A 


A 


0  ’ 

X 


9 

A 

W 

yi 

«/ 


3 


5 


than  the  rate  of  thermal  wave  penetration,  as  this  hot  layer  burns  out,  a  new 
thermal  layer  may  not  have  a  chance  to  be  established  for  sustained  burning. 
Thus,  propellant  sample  is  extinguished  after  the  first  thermal  layer  is 
burned  out.  The  chemical  kinetic  characteristics  for  solid  propellants  also 
have  significant  influence  on  the  momentary  ignition  phenomenon.  Since 
RDX/CTPB  propellant  requires  lower  activation  energy  than  HMX/CTPB,  it  is 
easier  to  establish  self-sustained  ignition  as  shown  in  Pictures  24  onwards. 
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III.  MECHANISM  OF  CRACK  PROPAGATION  IN  BURNING  SOLID  PROPELLANT 


3.1  Motivation  and  Objectives 

Increase  in  burning  surface  area  caused  by  propagating  cracks  can  affect 
the  actual  performance  of  a  rocket  motor  and  in  the  extreme  could  lead  to  de- 
f lagrat i on-to-detonat i on  transition  (DDT)  or  transition  to  detonation  due  to 
some  unresolved  mechanisms  (XDT).  The  study  of  propagating  cracks  in  solid 
propellants  under  closely  simulated  motor  operating  conditions  is  therefore 
important  in  understanding  the  off  design  operation  and  detonation  initiation 
phenomena  in  solid  propellants. 


One  of  the  major  objectives  of  this  research  program  is  to  gain  deeper 
understanding  of  the  mechanism  of  solid  propellant  grain  fracture  during 
rocket  motor  firing.  This  subject  is  of  special  interest  to  investigators  of 
rocket  motor  grain  integrity.  Numerous  studies  have  been  conducted  to  inves¬ 
tigate  crack  propagation  phenomena  by  employing  only  solid-mechanics  analyses 
and  experimental  procedures.  However,  very  little  effort  has  been  made  on 
crack  propagation  of  solid  propellants  under  burning  conditions.  It  is, 
therefore,  quite  obvious  that  the  study  of  crack  propagation  in  burning  solid 
propellants  would  be  most  realistic  and  instrumental  in  understanding  the  me¬ 
chanism  of  grain  fracture  due  to  crack  propagation  under  rocket  motor 
operating  conditions.  \ 

Specific  objectives  of  this  investigation  were: 


1.  To  develop  an  experimental  technique  to  study  crack  propagation  in  a 
composite  propellant  under  burning  conditions; 
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to  determine  crack  propagation  velocities  in  solid  propellant  sam¬ 
ples  under  burning  conditions; 

to  observe  any  abnormalities  of  the  propagating  crack; 

to  study  the  effects  of  chamber  pressure,  pressurization  rate,  and 
crack  sample  geometry  on  crack  propagation  velocity; 

to  analyze  the  stress  and  strain  state  at  the  crack  tip  by  utilizing 
a  finite  element  analysis; 

to  compare  the  observed  mass  burning  rate  with  that  calculated  from 
steady-state  burning  rate  law; 

to  modify  the  existing  test  rig  so  that  the  burning  can  be  interrup¬ 
ted  and  the  sample  can  be  recovered;  and 

to  postulate  and  identify  the  major  mechanism  of  dynamic  crack 
propagation. 

I  .  ‘ 

j 

3.2  Setup  of  Crack  Propagation  Test  Rigs 

Two  windowed  test  chambers,  one  small  and  one  large,  have  been  designed 
and  constructed  for  observing  the  crack  propagation  phenomena  in  burning 
solid  propellant  samples.  The  small  test  rig  was  adapted  from  the  chamber 
used  for  tip  ignition  studies.  The  schematic  diagram  and  a  part  of  the  re¬ 
sults  obtained  from  this  chamber  were  included  in  the  last  annual  report.  In 
order  to  study  the  crack  propagation  phenomena  under  more  suitable  condi¬ 
tions,  a  number  of  modifications  were  incorporated  into  the  design  of  a 
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larger  test  rig.  A  schematic  diagram  of  this  test  rig  is  shown  in  Fig.  10. 

Specific  features  and  improvements  of  the  second  test  rig  are: 

a)  increased  chamber  width; 

b)  increased  chamber  depth; 

c)  stronger  window  retainer  design; 

d)  more  complete  instrumentation  and  data  acquisition;  and 

e)  improved  safety  characteristics. 

With  the  increased  chamber  width,  the  sample  geometry  becomes  more 
flexible  and  the  influence  of  reflected  compression  waves  from  chamber  walls 
is  significantly  reduced.  Also,  the  crack  propagation  time  period  is  made 
long  enough  to  allow  more  detailed  observations  of  the  physical  process.  By 
increasing  the  chamber  depth,  possible  leaks  between  the  driving  motor  assem¬ 
bly  and  the  main  chamber  are  eliminated.  In  addition,  thicker  plexiglass 
windows  can  be  used  to  strengthen  the  window  assembly  so  that  the  possibility 
of  flame  spreading  between  the  propellant  sample  and  sacrificial  window  is 
substantially  reduced. 

A  photograph  of  the  main  chamber  and  window  retainer  is  shown  in  Fig. 
11.  The  window  assembly  consisting  of  a  1/2"  thick  sacrificial  window  and  a 
3"  thick  main  window  are  shown  together  with  the  window  retainer  in  Fig.  12. 

In  the  large  test  chamber,  two  window  openings  (a  circular  and  a  rectan¬ 
gular  one)  ware  made  in  the  window  retainer.  Ihe  purpose  of  designing  two 
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Fig.  10  Schematic  Diagram  of  Crack  Propagation  Chamber 
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viewing  openings  instead  of  a  single  long  one  is  to  avoid  window  buckling  and 
retainer  plate  warpage.  Through  the  circular  opening,  one  could  observe  the 
onset  of  the  combustion  product  gases  entering  the  test  chamber  from  the 
driving  motor. 

In  the  large  test  rig,  more  pressure  transducers  ware  used  for  recording 
p-t  traces  during  test  firings.  A  maximum  of  six  piezoelectric  transducers 
can  be  used  simultaneously.  One  of  these  transducers  could  be  mounted  on  the 
driving  motor  and  the  rest  mounted  on  the  main  chamber  at  the  crack  entrance, 
initial  crack  tip  region,  the  center  of  the  sample  beyond  the  initial  crack 
tip,  the  upper  cavity,  and  the  lower  cavity  (see  Fig.  10). 

Two  ports  for  mounting  break  wires  or  thermocouples  in  the  main  chamber 
were  designed  and  constructed  to  detect  the  arrival  of  the  hot  gases  pene¬ 
trating  the  propellant  sample  at  the  crack  tip.  Ports  were  made  in  the  free 
volune  regions  of  the  chamber  for  optical  illumination  of  the  sample  surfaces 
which  are  not  initially  in  contact  with  the  combustion  gases.  This  is  useful 
for  obtaining  information  of  the  instantaneous  boundaries  of  the  propellant 
sample  in  the  free  volume  region.  These  two  illumination  ports  can  also  be 
used  for  other  purposes  as  well;  for  example,  one  could  mount  transducers, 
photodiodes,  break  wires,  or  thermocouples  so  that  measurements  can  be  made 
in  the  free  volume  regions  or  the  upper  part  of  the  propellant  sample.  An 
additional  port  for  collecting  combustion  product  gas  samples  was  also  machi¬ 
ned  in  the  main  chamber  to  obtain  species  concentration  at  a  specified  time 
period  during  crack  propagation. 

An  assembly  of  safety  head,  consisting  of  a  specified  bursting  pressure 
diaphragm  and  a  retainer  is  installed  near  the  entrance  of  the  main  test 
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chamber  (see  Fig.  10).  This  setup  keeps  the  chamber  safe  in  case  of  over¬ 
pressurization.  It  is  also  advantageous  for  obtaining  higher  rates  of  de¬ 
pressurization,  since  this  can  help  the  recovery  of  the  propellant  sample  by 
quenching  the  burning  surfaces. 

In  sample  preparation,  a  propellant  block  is  first  machined  to  the  desi¬ 
red  dimensions.  It  is  then  clamped  in  between  two  plexiglass  tamplets  for 
trimming  into  a  final  configuration  (see  Fig.  13) .  A  3-mm  deep  slit  is  made 
by  a  sharp  razor  blade  at  the  crack  tip  to  help  the  initiation  of  the  crack 
growth  under  rapid  pressurization. 

It  is  important  to  note  that  shortly  after  the  flame  has  reached  the  tip 
of  the  narrow  slit,  the  subsequent  crack  propagation  beyond  the  original  tip 
of  the  slit  will  have  natural  crack  geometry  rather  than  the  artificial  geo¬ 
metry  of  the  initial  slit.  Also,  the  triangular  opening  of  the  propellant 
sample  at  the  entrance  is  to  facilitate  pressure  loading  on  the  slanted  sur¬ 
faces.  The  wide  triangular  opening  should  not  be  regarded  as  a  part  of  the 
initial  crack  tip  contour.  The  essential  part  of  crack  geometry  is  the  pre¬ 
fabricated  slit.  The  type  of  the  propellant  used  was  AP/HTPB  (73/27)  with  an 
average  AP  particle  size  of  200  m. 

Figure  14  shows  a  closeup  view  of  the  propellant  installed  in  the  test 
chamber.  The  sample  is  held  in  its  position  by  three  brass  retainers.  The 
plexiglass  window  assembly  is  designed  to  compress  the  sample  between  the 
rear  wall  of  the  chamber  and  the  window  to  avoid  any  leakage. 

A  photograph  of  the  assembled  test  rig  is  shown  in  Fig.  15.  In  this 
setup  the  driving  motor  is  mounted  below  the  horizontal  metal  plate  and 


Fig.  13  Photograph  of  the  Propellant  Sample  Clamped  Between 
Plexiglass  Tamplets 


cannot  be  seen  in  this  picture.  The  product  gases  from  the  driving  motor 
flow  upward  to  the  main  test  chamber  and  a  portion  of  the  hot  gases  enters 
the  crack  cavity. 

In  order  to  secure  the  test  rig,  several  angle  irons  are  used  to  bolt 
the  test  chamber  tightly  to  the  test  frame.  The  exhausted  gases  leave  the 
chamber  through  a  long  vertical  pipeline. 

3.3  Sequence  of  Events 

During  a  test  firing,  the  propellant  sample  cavity  is  pressurized  by 
high  temperature  and  high  pressure  gases  produced  from  the  small  driving 
motor.  The  above-mentioned  propellant  sample  configuration  permits  signifi¬ 
cant  displacements  of  the  crack  walls,  so  that  high  stresses  and  strains  can 
occur  at  the  tip  region  leading  to  crack  propagation.  An  electric  primer  is 
used  to  ignite  the  boost  propellant  and  the  igniter  charge.  The  combustion 
product  gases  from  the  driving  motor  flow  through  a  multi-perforated  nozzle 
into  the  main  chamber  and  enter  into  the  crack  cavity.  The  propellant  surfa¬ 
ces  in  the  cavity  become  ignited  shortly  after  the  penetration  of  the  hot 
gases.  The  gases  generated  from  the  burning  crack  surface  may  not  be  able  to 
leave  the  cavity  area  immediately,  therefore  causing  local  pressurization, 
mechanical  deformation,  and  crack  propagation  and  branching.  The  initial 
pressurization  of  the  crack  cavity  is  controlled  by  the  mass  and  energy 
fluxes  of  the  product  gases  coming  from  the  driving  motor.  The  pressuriza¬ 
tion  rates  generated  from  the  driving  motor  ranged  between  1  and  35  GPa/s, 
while  the  maximum  observed  pressure  in  the  test  chamber  was  80  MPa. 

3.4  Data  Acquisition  System 
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The  data  acquisition  system  consists  of  two  major  parts:  a)  a  pressure 
recording  system  and  b)  an  ignition  event  recording  system.  Pressure  measu¬ 
rements  are  taken  at  six  locations  described  previously.  The  pressure  signal 
from  a  transducer  is  sent  to  a  charge  amplifier  and  the  amplified  signal  is 
recorded  on  a  transient  waveform  recorder  and  a  high-speed  magnetic  tape  rec¬ 
order.  A  high-speed  16-mm  motion  picture  camera  is  used  to  film  the  ignition 
event  in  the  chamber.  The  maximum  framing  rate  of  the  camera  is  44,000  pic¬ 
tures  per  second.  The  detailed  discussion  of  the  data  acquisition  system  was 
presented  in  the  last  annual  report  and  hence  it  is  not  repeated  in  this 


report. 


srimental  Results 


In  order  to  study  the  effect  of  pressurization  rate  and  sample  geometry 
on  crack  propagation  processes,  other  variables  such  as  propellant  type,  ini¬ 
tial  temperature,  initial  angle  of  the  triangular  cavity,  were  held  constant 
in  this  investigation.  Three  sets  of  p-t  traces  obtained  in  test  firing  No. 
DNCP2-14,  DNCP2-15,  and  DNCP2-13  under  initial  pressurization  rates  of  10.4, 
8.3  and  33.4  GPa/s  are  shown  in  Figs.  16,  17,  and  18,  respectively. 

Figure  16  shows  the  set  of  p-t  traces  with  the  medium  pressurization 
rate  at  the  crack  tip.  These  traces  were  recorded  at  two  different  sampling 
rates  (2.5  ps  for  the  initial  uprising  portion  of  7.5  ms  and  2.5  ms  per 
sample  for  the  remaining  portion)  resulting  in  an  abrupt  change  in  pressure 
traces  near  pressure  spikes.  The  rapid  pressure  decay  after  the  peak  pres¬ 
sure  is  due  to  the  rupture  of  the  bursting  diaphragm.  As  can  be  seen  from 
Fig.  16,  the  first  discernible  pressure  rise  at  the  sample  center  G4  (which 
is  beyond  the  initial  crack  tip)  occurs  later  than  that  at  the  crack  tip  G3. 
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Fig.  16  experimental  data  of  crack  propagation  test  firing  no,  DNCP2-14 

(framing  rate  =  16,000  pps) 


Fig.  17  experimental  data  of  crack  propagation  test  firing  no.  DNCP2-15 

(framing  rate  =  16,000  pps) 


Fig.  18  experimental  data  of  crack  propagation  test  firing  no.  DNCP2-13 

(framing  rate  =  32,000  pps) 
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However,  the  pressure  and  the  pressurization  rate  of  <34  is  higher  than  those 
of  G3.  This  higher  pressure  and  pressurization  rate  at  <34  is  believed  to  be 
caused  by  large  specific  surface  area  generated  by  microstructure  damage  of 
the  propellant  during  crack  propagation  and  branching  as  seen  from  the  high¬ 
speed  film  to  be  discussed  in  a  later  section. 

The  pressure-time  traces  in  Fig.  12,  corresponding  to  the  lower  initial 
pressurization  rate,  show  that  the  maximum  pressure  near  the  center  of  the 
sample  is  lower  than  that  at  the  tip.  After  the  initial  time  period  of  3  ms, 
the  pressurization  rates  at  (34,  G3,  and  G2  become  almost  the  same.  Comparing 
Fig.  17  (low  initial  pressurization  rate)  with  Fig.  16  (medium  initial  pres¬ 
surization  rate)  it  can  be  seen  that  the  pressurization  rate  is  increasing  in 
the  downstream  direction  for  the  experiment  with  the  higher  initial  pressuri¬ 
zation  rate. 

Figure  18  reveals  an  important  phenomenon  occuring  in  the  firing  test 
with  the  highest  initial  pressurization  rate  at  <33.  The  pressure  at  the 
center  of  the  sample  rose  so  abruptly  and  became  saturated  in  a  very  short 
time.  From  this  figure,  one  can  expect  that  the  center  of  the  sample  must  be 
burning  severely  in  order  to  have  such  high  gasification  and  pressurization 
rates.  Comparing  p-t  traces  of  Fig.  18  with~those  of  Figs.  16  and  17,  it  can 
be  observed  that  the  response  of  <34  is  steeper  for  higher  initial  pressuriza¬ 
tion  rate  at  G3.  One  explanation  for  this  observation  is  the  increased 
degree  of  damage  near  the  crack  tip  region  due  to  higher  pressure  loading. 
This  damage  can  contribute  significantly  to  provide  a  larger  burning  surface 
area  and  consequently  increase  the  local  mass  burning  rate. 
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In  tnese  tests,  the  fuse-wire  signals  shown  on  Figs.  16,  17,  and  18  were 
not  helpful  in  giving  information  for  the  arrival  of  crack  front.  This  is 
due  to  the  fact  that  fuse  wires  were  mounted  too  close  to  the  tip  where  large 
deformation  occurred.  New  locations  for  fuse  wires  in  the  upper  part  of  the 
propellant  sample  are  planned  for  future  firing  tests. 

The  crack  propagation  event  was  recorded  by  means  of  high-speed  photog¬ 
raphy.  Two  typical  film  records  corresponding  to  test  firings  DNCP2-14  and 
DNCP2-15  are  shown  in  Figs.  19  and  20,  respectively.  As  one  can  see  from 
these  films,  hot  gases  from  the  driving  motor  illuminated  the  crack  cavity 
and  outlined  the  crack  walls.  In  these  figures,  each  vertical,  column  of 
light  is  the  outlined  crack  at  a  different  time.  The  event  proceeds  from 
left  to  right,  with  a  time  interval  of  59  ys  between  consecutive  pictures. 
The  event  shown  in  Fig.  19  is  described  as  follows. 

Observations  and  Interpretation 
Hot  igniter  gas  starts  to  enter  the  triangular 
cavity. 

A  typical  mode  I  crack  propagation  occurs  with 
the  crack  tip  displacement  in  the  axial 
direction.  Near  the  crack  tip  a  small 
triangle  is  formed  indicating  the  existence 
of  only  one  crack.  The  dark  horizontal  line 
is  due  to  the  view  blockage  of  the  window 
retainer  material  between  the  two  viewing 
ports. 

A  dramatic  change  takes  place  in  the 
propagation  mechanism.  Two  spherical  regions 
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are  developing  near  the  tip  and  a  significant 
enlargement  of  the  luminous  volume  occurs. 

This  shows  that  branching  starts  at  two 
different  places  located  near  the  shoulders 
of  the  luminous  volume  and  propagates  outward. 

8-10  413-531  A  mushroom-shaped  luminous  region  is 

developing  near  the  tip  of  the  crack.  In  the 
meantime,  the  crack  walls  became  almost 
parallel,  and  the  crack  tip  itself  loses  its 
definite  triangular  shape  and  became  a  jagged 
horizontal  line.  One  explanation  for  the 
mushroom-shaped  region  and  the  corresponding 
change  in  crack  geometry  is  the  microstructure 
damage  which  permits  the  penetration  of  hot 
gases  into  the  damaged  region  near  the  crack 
tip. 

11  590  A  reduction  of  the  mushroom-shaped  region  can 

be  observed.  The  severe  burning  near  the 
burning  near  the  crack  tip  can  cause 
mechanical  compression  of  the  material  ahead 
of  the  jagged  horizontal  front.  This 
compression  may  lead  to  the  close  up  of  some 
micropores,  and  hence  reduce  the  luminosity 
of  the  mushroom-shaped  region. 


12-13  649-708  The  continued  burning  and  gas  penetration  at 

the  tip  front  enlarges  the  luminous  zone. 

14  767  Same  pores  of  the  burning  propellant  are 

extended  far  enough  to  allow  the  penetration 
of  the  flame  to  the  free  volume  regions, 
causing  illumination  of  the  outer  boundary 
surfaces  of  the  sample.  Large  mechanical 
deformation  of  the  outer  boundaries  is 
observed. 

15-21  826-1180  Severe  burning  is  taking  place  near  the  center 

of  the  test  sample.  The  mechanical 
deformation  becames  so  large  that  the  free 
volumes  almost  disappear.  The  mushroom-shaped 
region  continues  to  enlarge  at  the  center  of 
the  sample.  The  total  mass  burning  rate 
increases  significantly  with  time. 

Figure  20  shows  the  crack  propagation  event  under  a  lower  initial  pres¬ 
surization  rate  at  G3.  The  general  sequence  of  the  event  is  similar  to  that 
of  Fig.  19.  The  major  differences  between  the  two  events  are: 

1)  The  bifercation  of  the  spherical  regions  starts  much  later  for  the 
lower  pressurization  case.  The  spherical  regions  are  not  as  lumi¬ 
nous  as  those  shown  in  Fig.  19. 


60 


2)  The  development  of  a  mushroom-shaped  region  also  occurs  at  a  later 
time. 

3)  The  flame  penetration  through  the  sample  into  the  free  volume  re¬ 
gions  also  occurs  at  a  later  time. 

4)  Opposed  to  Fig.  19,  unsynmetric  mushroom-shaped  regions  are 
observed.  This  is  probably  caused  by  the  close  up  of  some  micropo¬ 
res  due  to  uneven  mechanical  deformation. 

The  relationship  between  vertical  crack-tip  displacement  and  chamber 
pressure  was  achieved  by  combining  the  pressure-time  traces:  and  the  dis¬ 
placement  vs.  time  observations.  Examination  of  this  relationship 

j 

illustrates  the  influence  of  the  pressurization  rate  and  crack  sample  geo- 

1 

metry  on  crack  propagation  yelocity.  Displacement  vs.  pressure  curves  for 

i 

three  tests  having  the  sam4  initial  sample  geometry,  but  different  pressuri- 

i 

ration  rates,  are  shown  in  fig.  21a.  A  linear  approximation  was  used  for  the 
data  points  obtained  in  each  test.  It  is  obvious  from  this  figure  that  the 
tests  resulted  in  three  linear  relationships  between  displacement  and  pres- 

I 

sure,  and  that  the  slopes  of  the  displacement  vs.  pressure  curves  decrease 
with  increasing  pressurization  rate.  Considering  that  only  pressurization 
rate  and  sample  geometry  were  varied,  the  relationship  between  the  slopes  of 
the  displacement-pressure  curves  and  these  parameters  can  be  generally  writ¬ 
ten  as: 

(31) 


ds 

dp 


F(^  >  geometry) 


PRESSURE,  p,  MPa 


Fig.  21a  Measured  Crack  Tip  Displacement  vs.  Chamber  Pressure  for 
Crack  Samples  with  the  Same  Initial  Geometries 
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For  the  tests  shown  in  Fig.  21a,  where  the  geometry  was  held  constant, 
the  above  expression  reduces  to  the  following: 

~  -  f  (32) 

ap  at 

This  function  was  approximated  by  a  power  law,  whose  coefficients  were  deter¬ 
mined  by  a  linear  approximation  of  the  log-log  plot  of  ds/dp  versus  dp/dt  for 
the  tests  with  the  fixed  geometry.  The  function  was  found  to  be  closely  ap¬ 
proximated  by  the  following  relationship 

-  0.086  ($)"  0,363  (33) 

dp  dt 

The  above  equation  was  used  to  derive  the  following  experimental  correlation 
relating  crack  propagation  velocity  and  pressurization  rate: 

0.637 

(34) 

It  can  be  seen  that  the  derived  correlation  is  in  agreement  with  the  afore¬ 
mentioned  observation  that  the  crack  propagation  velocity  increases  as  the 
chamber  pressurization  rate  becomes  higher. 

To  examine  the  influence  of  different  propellant  sample  geometries  on 
observed  crack  propagatoin  velocities,  three  crack-tip  displacement  vs.  pres¬ 
sure  curves  are  plotted  in  Fig.  21b.  During  the  three  separate  tests  with 
different  sample  geometry,  the  pressurization  rate  was  held  constant,  ap¬ 
proximately  15,000  MPa/s.  The  sample  geometry  differed  in  the  taper  angle  of 
the  external  surfaces  of  the  propellant  sample,  which  changes  the  free  volume 
regions  of  the  test  chamber.  The  sample  with  the  maximan  free  volume  corres¬ 
ponds  to  the  curve  with  the  highest  ds/dp  shown  in  Fig.  21b.  Considering  the 
fact  that  these  tests  were  conducted  under  the  same  average  pressurization 


^  ■  Vc  (m/s)  «  0.086  (MPa/s) 


DISPLACEMENT  OF  LUMINOUS  CRACK-TIP  FRONT, 
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rate,  it  can  be  easily  found  that  the  steeper  the  slope  of  the  curve  in  Fig. 
21b,  the  higher  the  average  crack  propagation  velocity. 

Combined  data  obtained  from  several  tests  with  different  geometries  and 
pressurization  rates  are  shown  in  Fig.  22.  The  solid  line  is  the  correlation 
given  by  Eg.  (34) ,  and  the  dashed  lines  represent  the  similar  relationships 
at  different  taper  angles.  Since  the  data  points  of  these  curves  are  limit¬ 
ed,  the  main  usefulness  of  Fig.  22  is  to  display  the  qualitative  trends  of 
the  effect  of  dp/dt  and  taper  angle  on  crack  propagation  velocity. 

In  order  to  explain  the  crack  propagation  velocity  dependence  on  pres¬ 
surization  rate  and  sample  geometry,  two  cornnon  propagation  mechanisms  were 
examined.  The  first  mechanism  considers  that  the  crack  is  propagating  at  its 

terminal  velocity,  Vfc,  which  is  approximately  equal  to  one-third  of  the  speed 

33 

of  sound  in  the  propellant  as  suggested  by  Swanson .  As  the  pressurization 
rate  is  increasing,  the  dynamic  and  relaxation  moduli  of  the  viscoelastic 
propellant  are  also  increasing.  This  results  in  a  higher  terminal  velocity 
in  the  propellant,  according  to  the  following  formula: 

Vt  -  i  /e/ps  (35) 

where  E  is  the  dynamic  modulus  and  pg  is  the  density  of  the  propellant 

sample.  Over  the  range  of  pressurization  rates  tested,  johanges  in  moduli  are 

approximately  10%;  therefore,  these  small  changes  in  moduli  do  not  result  in 

significant  increase  in  crack  propagation  velocity  as  measured  in  this  work. 

34  1 

This  observation  has  been  made  by  Gent  end  Marteny  in  their  work  with  natu¬ 
ral  rubbers  under  rapid  loading  rates.  The  first  mechanism  does  not  consider 
any  changes  in  the  crack  propagation  velocity  due  to  different  sample  geome- 
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Fig.  22  The  Combined  Effect  of  Pressurization  Rate  and  Sample 
Geometry  on  Crack  Propagation  Velocity 
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tries. 

The  second  crack  propagation  mechanism  to  be  discussed  considers  that 
the  velocity  of  the  running  crack  is  limited  by  the  speed  of  the  deforming 
crack  contour.  In  order  to  have  crack  propagation,  a  sufficient  stress  state 
must  be  maintained  at  the  crack  tip.  For  geometries  tested  in  this  study, 
the  large  deformation  of  the  crack  cavity  observed  from  the  high-speed  films 
is  believed  to  provide  the  high  stress  distribution  near  the  crack  tip.  The 
crack  tip  front  shown  by  the  luminous  contour  seems  to  be  limited  by  the  rate 
of  deformation  of  the  crack  cavity  walls.  However,  actual  crack  tip  propag¬ 
ation  velocity  may  exceed  the  observed  velocity  obtained  from  the 
displacement  cf  the  luminous  front.  The  microstructural  damage  of  the  pro¬ 
pellant  ahead  of  the  luminous  front  may  not  be  seen  from  high-speed  films. 
Sane  of  the  micro  cracks  are  submerged  beneath  the  surface,  and  also  there  is 
a  time  lag  for  the  gas  to  penetrate  these  micro  cracks. 


The  quasistatic  application  of  this  mechanism  would  predict  that  the 
crack  propagation  velocity  is  proportional  to  the  first  power  of  pressuriza¬ 
tion  rate.  This  is  because  of  the  fact  that  a  given  static  pressure  load 
would  correspond  to  a  certain  displacement  of  the  crack  contour.  This  should 
result  in  the  three  curves  in  Fig.  21a  coalescing  into  one.  However,  the  ex¬ 
perimental  results  of  Fig.  21a  show  a  significant  dependence  of  ds/dp  on 
dp/dt.  The  changes  of  slope  at  different  pressurization  rates  shown  in  this 
figure  sure  believed  to  be  caused  by  the  delay  of  the  sample  leg  displaanent 
at  higher  loading  rates  due  to  the  effect  of  inertia  during  the  acceleration 
of  the  sample  mass. 
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From  the  above,  one  can  deduce  that  the  dynamic  application  of  this  me¬ 
chanism  would  show  weaker  dependence  of  V  on  dp/dt  than  the  quasistatic  one. 
The  exponential  coefficient  of  dp/dt  in  the  dynamic  application  would  be 
lower  than  one  which  is  in  agreement  with  the  experimental  results  of  Eq. 
(34) . 

Similarly,  the  changes  in  slope  at  different  geometries  shown  in  Fig.  22 
are  believed  to  be  caused  by  the  decrease  in  resistance  of  the  sample  leg 
displacement  for  larger  taper  angles.  This  implies  that  the  greater  the  free 
volume,  the  quicker  the  response  of  the  propellant  sample  to  the  pressure 
loading  and  consequently  the  higher  the  crack  propagation  velocities.  The 
approximate  magnitudes  and  the  tendencies  of  the  curves  in  Fig.  21a  and  21b 
were  simulated  by  solving  analytically  a  very  simplified  model  consisting  of 
a  Voigt-Kelvin  parallel  spring-damper  element  in  series  with  a  point  mass. 
The  coefficients  for  the  system  were  approximated  from  results  of  the  static 
finite  element  analysis  and  sample  geometry.  The  results  of  the  above  simple 
analysis  produce  qualitatively  the  same  changes  in  ds/dp,  due  to  the  differ¬ 
ent  loading  rates  and  geometries. 

It  is  believed  that  at  much  higher  pressurization  rates  than  those 
tested  in  this  study,  higher  crack  propagation  velocity  would  be  reached. 
Under  these  conditions,  the  influence  of  sample  geometry  would  be  less  signi¬ 
ficant  and  the  first  crack  propagation  mechanism  would  be  important. 

33  34 

Swanson  and  Gent  and  Marteny  measured  crack  propagation  velocities 
in  inert  propellant  (Solithane  113)  and  rubber,  respectively.  Swanson  repor¬ 
ted  the  highest  observed  velocities  of  approximately  40  Vs,  while  Gent  and 
34 

Marteny  measured  highest  velocities  of  80  Vs  in  specimens  with  no  pre- 
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imposed  strains  in  the  crack  direction.  It  is  interesting  to  note  that  their 
measured  highest  velocity  is  in  the  same  order  of  magnitude  as  the  observed 
1 'mi nous  front  propagation  velocity  in  burning  solid  propellants.  The  actual 
crack  propagation  velocity  can  be  substantially  higher  than  that  of  the  lumi¬ 
nous  front  for  reasons  stated  above. 

3.6  Finite  Element  Analysis 

An  analytical  stress  analysis  of  the  propellant  sample  was  performed 
using  a  finite  element  techique.  The  objectives  of  this  analysis  were: 
a)  to  understand  how  the  sample  would  deform  in  a  cold  environment  under 
uniform  internal  pressure  loading;  b)  to  predict  the  way  in  which  the  pro¬ 
pellant  sample  fractures;  and  c)  to  estimate  the  contribution  of  mechanical 
deformation  in  the  overall  measured  crack  tip  displacement. 

The  propellant  sample  geometry  modeled  was  similar  to  that  described  in 
Section  3.2.  A  finite  radius  of  1.78  mm  was  selected  at  the  crack  tip  in 
order  to  use  a  manageable  number  of  mesh  elements  to  study  the  local  stress 
and  strain  distributions.  The  finite  radius  approximation  reduces  the  accu¬ 
racy  of  the  solution  near  the  tip.  However/  the  solution  is  useful  to  have  a 
qualitative  representation  of  the  mechanical  deformation  of  the  propellant 
sample  near  the  crack  tip. 

The  employed  finite  element  grid  was  generated  by  using  an  interactive 
mesh  generation  program  and  the  grid  pattern  was  shown  in  the  last  annual 
report.  The  material  properties  assumed  for  the  model  where  an  elastic  modu¬ 
lus,  E,  equal  to  55  MPa  and  Poisson's  ratio,  v,  equal  to  0.5.  Even  though 
elastic  behavior  was  considered,  the  viscoelastic  character  of  the  sample  was 
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approximated  by  substituting  the  relaxation  modulus  for  the  elastic  modulus. 
According  to  the  correspondence  principle,  the  selection  of  the  relaxation 
modulus  was  based  upon  the  test  pressurization  rate  and  initial  temperature 
of  the  propellant.  This  is  a  commonly  used  procedure  in  solid  propellant 
grain  deformation  analysis. 

A  static  analysis  of  the  described  model  was  solved  with  a  nonlinear 
finite  element  code  (ABAQUS).^  The  program  permitted  incompressible  mater¬ 
ial  behavior  and  large  deformations  which  were  needed  in  order  to  solve  the 
problem.  The  obtained  results  are  discussed  in  terms  of  the  following  three 
subject  areas  of  interest:  a)  displacement  of  the  crack  tip  along  the  crack 
centerline;  b)  contour  deformation;  and  c)  stress  and  strain  distributipn 
near  the  crack  tip. 

a)  Crack  Tip  Displacement 

In  terms  of  crack  tip  displacement,  the  pressure  loading  on  the  inner 
crack  walls  resulted  in  a  mechanical  displacement  along  the  crack  centerline. 
The  calculated  vertical  displacement  was  due  to  the  fact  that  the  free  boun¬ 
daries  of  the  sample  allowed  for  horizontal  expansion.  The  displacements  of 
the  crack  tip  node  were  determined  for  several  pressurization  rates,  and  the 
calculated  results  indicated  the  nonlinearity  of  the  problem.  The  results  of 
mechanical  displacement  vsl.  pressure  can  be  easily  converted  to  displacement 
vs.  time  data,  assuming  that  the  pressure-time  relationship  is  given  and  the 
quasisteady  assumption  is  Valid. 

TVfo  mechanical  displacement  curves  obtained  in  this  way  are  shown  in 
Fig.  23;  one  corresponding  to  a  low  pressurization  rate  of  3,000  MPa/s  and 
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the  other  to  a  high  pressurization  rate  of  30,000  MPa/s.  In  order  to  deter¬ 
mine  the  contribution  of  the  mechanical  deformation  to  the  overall  measured 
luminous  front  displacement  at  the  crack  tip,  the  experimentally  observed 
displacesnents  under  these  two  pressurization  rates  were  plotted  in  the  same 
figure.  In  the  low  dp/dt  case,  it  can  be  seen  that  the  mechanical  dis¬ 
placement  was  considerably  snaller  than  the  overall  displacement.  On  the 
other  hand,  in  the  high  dp/dr  case,  a  more  significant  contribution  of  the 
mechanical  deformation  the  overall  luminous  front  displacement  at  the 
crack  tip  can  be  observed.  The  conclusion  is  that  the  higher  the  pressuriza¬ 
tion  rate,  the  larger  the  contribution  of  the  mechanical  deformation  to  the 
overall  luninous  front  displacement  at  the  crack  tip.  In  order  to  determine 
the  accuracy  of  the  above  results,  it  is  necessary  to  examine  the  quasisteady 
assumption  made  in  the  finite  element  analysis.  A  simple  lutped-mass  analy¬ 
sis  demonstrates  that  time  lags  due  to  inertia  effects  are  of  the  same  order 
of  magnitude  as  the  duration  of  the  experiment.  Therefore,  the  higher  the 
leading  rate,  the  larger  the  errors  due  to  the  quasisteady  assumption.  In 
addition,  the  static  solution  gives  an  upper  boundary  of  the  displacement  at 
a  given  time,  since  the  inertia  effect  causes  actual  displacement  to  lag 
behind  the  calculated  value. 

b)  Contour  Deformation 

Wnen  fracture  at  the  crack  tip  is  not  allowed,  the  mechanical  deform¬ 
ation  at  the  crack  tip  tends  to  slightly  round  the  original  triangular  shape 
of  the  crack  inner  wall.  However,  the  calculated  crack  tip  deformation  is 
snail,  as  shown  in  Fig.  23,  when  compared  with  the  observed  high  rates  of 
contour  changes  from  experiments.  This  implies  that  the  crack  propagation 
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must  be  caused  by  fracture.  Finite  element  solution  showed  that  the  outer 
surface  of  the  sample  deforms  towards  the  chamber  wall  as  anticipated.  How¬ 
ever,  if  no  crack  propagation  is  allowed,  the  finite  element  solution  will 
show  a  small  amount  of  deformation  of  the  propellant  sample  legs  toward  the 
chamber  wall.  If  one  or  more  cracks  were  propagating  the  leg  would  be  free 
to  deform  out  towards  the  chamber  wall,  creating  abrupt  changes  in  the  crack 
cavity  contour, 

c)  Str*:  s  and  Strain  States  Near  the  Crack  Tip 

In  ordfc.  to  understand  the  mode  of  failure  of  the  propellant  sample,  a 
typical  propellant  grain  failure  criterion  should  be  used.  An  experimentally 
derived  failure  envelope,  based  upon  the  maximum  principal  deviatoric  stress 
and  the  corresponding  strain,  is  a  cannon  criterion  for  determining  a 
loading-rate  independent  failure  of  solid  propellant  grains  subject  to  inter¬ 
nal  pressurization.  According  to  this  criterion,  the  propellant  will  begin 
to  crack  at  the  surface  near  a  point  where  the  stress  and  strain  states 
exceed  the  boundary  of  the  envelope.  The  two  important  characteristics  of 
failure  criterion  are:  1)  failure  is  dependent  upon  deviatoric  stress  and 
not  upon  hoop  or  maximum  principal  stress;  and  2)  the  criterion  is  based 
upon  the  level  of  the  combined  state  of  stress  and  strain. 

'*he  deviatoric  stress  contours  near  the  crack  tip  region  are  shown  in 
Fig.  24a.  This  figure  clearly  shows  that  the  maximum  deviatoric  stress  con¬ 
centration  is  near  the  shoulder  of  the  crack  tip  radius.  Four  locations  on 
the  crack  radius  are  indicated  on  the  figure  as  points  A,  B,  C,  and  D.  A 
plot  of  the  maximum  deviatoric  stress  vs.  the  corresponding  strains  for  these 
four  points  is  shown  in  Fig.  24b.  It  can  be  seen  that  the  maximum  deviatoric 
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Fig.  24a  Isopleths  for  Deviatoric  Stresses 


MAXIMUM  DEVIATORIC  STRESS,  <r,  MPa 
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Fig.  24b  Maximum  Deviatoric  Stress  vs  Corresponding  Strain  for  Crack 
Tip  Elements  at  a  Pressure  of  5.17  MPa 
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stress  occurs  at  D,  on  the  shoulder  of  the  crack  tip,  not  at  the  tip  zenith. 
In  addition,  the  maximum  strain  occurs  at  C.  This  implies  that,  for  the  geo¬ 
metry  and  load  conditions  modeled  by  the  finite  element  analysis,  the  most 
likely  area  for  failure  would  be  in.  the  vicinity  of  C  and  D,  and  not  at  the 
crack  tip.  The  symmetry  of  the  model  implies  that  the  failure  at  the  shoul¬ 
der  could  create  two  running  cracks.  In  this  case,  the  two  cracks  would  be 
extremely  close  together  at  the  onset  of  failure,  because  of  the  small  radius 
of  the  crack  tip.  The  consideration  of  dual  fracture  helps  to  explain  the 
transition  to  the  two  spherical  regions  and  then  to  a  square-shaped  geometry 
as  observed  on  the  film  during  high  pressurization  rate  tests.  A  multiple 
fracture  of  this  type  would  cause  dramatic  displacement  of  the  crack  cavity 
walls  toward  the  chamber  wall. 

3.7  Apparent  Versus  Steady-State  Mass  Burning  Rates 

An  estimation  of  the  apparent  mass  burning  rate  during  the  crack  pro- 
pagaton  can  give  valuable  information  about  the  degree  of  fracture  due  to 
crack  propagation  and  branching,  since  direct  measurement  of  the  burning  rate 
under  these  conditions  is  very  difficult.  One  method  of  obtaining  the  appar¬ 
ent  burning  rate  is  by  using  the  average  pressurization  rate  and  initial 
crack  geometry  as  given  in  the  following.  According  to  mass  continuity  for 
the  gases  in  the  chamber,  we  have 

rl  M  •  •  • 

=  ni  +  ra.  -  m  (36) 

dt  b  in  out 

After  a  short  initial  pressurization  period,  ir^  becomes  small  anc  can  be  ne¬ 
glected.  The  flow  at  the  exit  nozzle  can  be  considered  choked  because  of 
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high  chamber  pressure.  Using  these  assumptions  together  with  the  perfect  gas 
law,  Eq.  (36)  becomes: 


"b 


/R  T 


RT  dt 


c 


(37) 


where  dP^dt  represents  the  average  chamber  pressurization  rate.  In  the 
above  simplified  analysis,  both  Tc  and  were  treated  as  known  constants. 

Figure  25  shows  the  ratio  of  the  computed  apparent  mass  burning  rate 
(n^)  to  the  steady-state  mass  burning  rate  (m^  gg)  evaluated  at  the  average 
chamber  pressure  during  crack  propagation.  It  is  useful  to  note  that  the  ap¬ 
parent  mass  burning  rates  are  about  two  orders  of  magnitude  higher  than  the 
steady-state  mass  burning  rate.  The  ratio  increases  for  higher  pressuriza¬ 
tion  rates.  Also,  this  ratio  becomes  higher  for  test  chamber  having  larger 
free  volume  for  sample  deformation. 

This  apparent  increase  in  the  mass  burning  rate  is  believed  to  be  mainly 
caused  by  the  propellant  fracture  c.nd  microstructure  damage.  This  is  also 
supported  by  the  luminous  mushroom-shaped  region  observed  above  the  crack 
tip.  The  dynamic  (transient)  burning  effect^®  would  also  contribute  to  the 
total  mass  burning  rate.  The  relative  contribution  of  these  two  effects  is 
difficult  to  determine.  This  greatly  enhanced  mass  generation  rate  in  highly 
confined  regions  could  contribute  significantly  to  the  onset  of  detonation. 


RATIO  OF  MASS  BURNING  RATES,  *»/*».  •• 


Fig.  25  Effect  of  Pressurization  Rate  on 
the  Ratio  of  Apparent  to  Steady- 
State  Mass  Burning  Rate  during 
Crack  Propagation 
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3.8  Designed  Modifications  for  Achieving  Interrupted  Burning  During 

Crack  Propagation 

A  test  rig  modification  is  being  performed  in  order  to  recover  the  pro¬ 
pellant  crack  sample  from  interrupted  crack  propagation  experiments.  The 
propellant  sample  recovery  is  important  in  obtaining  direct  evidence  of  the 
microstructure  damage  and  in  achieving  better  understanding  of  the  crack  pro¬ 
pagation  and  branching  mechanisms  under  burning  conditions.  Careful 
considerations  were  given  to  the  test  rig  modification,  so  that  the  sample 
can  be  extinguished  at  a  prescribed  pressure  level  by  rapid  depressurization 
followed  with  inert  gas  injection. 

A  schematic  diagram  of  the  modified  test  rig  and  the  control  circuit  for 
depressurization  and  inert  gas  injection  is  shown  in  Fig.  26.  The  procedure 
to  achieve  the  propellant  sample  quenching  is  as  follows: 

1.  As  soon  as  the  pressure  in  the  chamber  reaches  a  prescribed  value  p^ 
(shewn  in  Fig.  27),  a  trigger  device  will  activate  the  primer  located  in 
the  second  igniter  assembly  at  time  t^. 

2.  The  product  gases  from  the  second  igniter  assembly  (see  Fig.  25)  will 
pressurize  the  free  space  between  the  two  bursting  diaphragms  and  rupture 
them  at  time  t^ *  to  cause  rapid  depressurization  of  the  test  chamber. 

3.  Immediate  injection  of  inert  gases  at  time  t2'  upon  a  time  delay  of  (t2'- 
t2) ,  where  t2  is  the  time  when  the  inert  gas  inlet  valve  was  activated. 
The  time  distance  t  between  t2  and  t^  can  be  preselected  and  adjusted  on 


a  timer 
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Fig.  26  Modified  Test  Rig  for  Sample  Recovery 


4.  The  inert  gas  injected  into  the  chamber  will  quench  the  inner  surfaces  of 
the  test  chamber  and  further  interrupt  the  burning  of  the  propellant 
sample. 

The  recovered  propellant  sample  will  be  examined  in  detail  under  a 
microscope  to  determine  the  degree  of  damage  encountered  during  the  burning. 
The  recovered  sample  will  be  extremely  useful  for  understanding  the  crack 
propagation  and  branching  mechanism.  One  preliminary  attempt  was  made 
towards  this  direction  without  using  the  above-mentioned  circuit;  and  a  small 
portion  of  the  propellant  was  recovered.  It  is  believed  that  larger  recover¬ 
ed  samples  will  be  obtained  after  the  test  rig  modification  is  completed. 
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IV.  SUMMARY  AND  CONCLUSIONS 

Some  important  observations  and  conclusions  from  the  tip  ignition  study 
are  summarized  as  follows: 

(1)  Ignition  of  nitramine  composite  propellants  was  studied  at  rapid 
pressurization  environment.  For  each  test  firing,  two  different  pro¬ 
pellant  samples  can  be  tested  simultaneously  under  the  same  flow 
conditions.  Three  independent  data  acquisition  systems  were  used  for 
detecting  the  onset  of  gas  evolution.  These  systems  include  the  ex¬ 
ternal  photodiode  system,  near-infrared  photodiode  system,  and 
high-speed  filming  system.  Consistent  results  were  obtained  from 
these  systems  in  determining  the  onset  of  gas  evolution. 

(2)  A  theoretical  model  with  chemical  kinetic  is  proposed  for  the  igni¬ 
tion  of  nitramine  composite  propellants.  These  kinetic  mechanisms 
and  data  are  helpful  in  the  interpretation  of  experimental  results. 

(3)  The  time  required  for  the  onset  of  gas  evolution  of  nitramine  compo¬ 
site  propellants  decreases  drastically  as  the  pressurization  rate  is 
increased. 

(4)  Nitramine  composite  propellants  are  much  me  re  difficult  to  ignite  as 
compared  to  the  conventional  AP-based  composite  propellants. 

(5)  With  the  same  binder  HMX-based  composite  propellants  were  found  to 
be  more  difficult  to  gasify  than  PDX-based  composite  propellants. 
This  is  probably  due  to  the  fact  that  both  the  thermal  diffusivity 
and  activation  energy  for  decomposition  of  RDX  are  lower  than  those 
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of  HMX-based  composite  propellants. 

(6)  Two  HMX-based  composite  propellants  with  different  binders  (CTPB  and 
HTPB)  did  not  show  significant  differences  in  the  time  required  for 
the  onset  of  gas  evolution.  This  is  probably  due  to  the  fact  that 
the  thermal  diffusivity  effect  and  the  chemical  kinetic  effect  of  de¬ 
composition  of  these  two  binders  tend  to  cancel  each  other. 

Several  important  observations  and  conclusions  obtained  from  the  crack 
propagation  studies  are  sixnmarized  below. 

(1)  Luminous  front  propagtaion  velocities  of  20  to  80  m/s  were  observed 
in  a  burning  propellant  crack  sample. 

(2)  For  high  pressurization  rates,  the  luminous  crack  tip  contour  under¬ 
goes  a  highly  discernible  transition  process  that  a  triangular 
contour  develops  into  a  mushroom-shaped  region.  The  evolution  of 
the  mushroom-shaped  region  is  believed  to  be  caused  by  gas  penetr¬ 
ation  and  flame  propagation  into  micropores  resulting  from  crack 
branching. 

(3)  An  experimental  correlation,  showing  strong  dependence  of  crack  pro¬ 
pagation  velocity  on  pressurization  rate  and  crack  sample  geometry, 
has  been  developed. 

(4)  The  apparent  mass  burning  rate,  deduced  from  the  measured  p-t 
traces,  was  found  to  be  approximately  two  orders  of  magnitude  higher 
than  the  steady-state  mass  burning  rate.  This  implies  that  the  pro¬ 
pellant  structure  damage  due  to  crack  propagation  and  branching  can 
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significantly  increase  the  burning  surface  area.  This  may  cause  the 
onset  of  detonation  in  strongly  confined  environments. 
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